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Rain and ice crystals depolarize radio waves along earth-satellite 
propagation paths. The magnitude of this depolarization is a function 
of incident polarization angle and is minimized when polarization and 
depolarizer symmetry axes coincide. A technique is presented for a 
direct determination of the medium's attenuation and depolarization 
for any incident polarization, based on measurements taken at two 
orthogonal polarizations. Some sample results from this technique are 
presented, using data collected at Crawford Hill, New Jersey using the 
19-GHz COMSTAR satellite beacon. 

I. INTRODUCTION 

Depolarization caused by rain and ice crystals is an important factor 
in the design of future satellite communication systems operating at 
frequencies above 10 GHz. 12 These systems will likely use dual, or- 
thogonal polarizations to increase transmission capacity. Quantitative 
knowledge of depolarization is necessary for determining the pair of 
polarizations that experience the least depolarization, for determining 
whether the isolation between any two polarizations is adequate during 
rain and ice depolarizing conditions, or for guiding the design of circuits 
for canceling crosstalk resulting from depolarization if the isolation is 
not adequate. 

Only a few measurements have been made of depolarization along 
earth-space propagation paths. 3-5 These measurements were not nec- 
essarily made at incident polarizations that produce minimum depo- 
larization. Aerodynamic forces acting on vertically falling raindrops are 
expected to orient raindrop symmetry axes vertical and horizontal on 
the average. Thus, rain depolarization is expected to minimize at linear 
horizontal and vertical polarizations. Maximum depolarization is ex- 
pected for linear polarizations oriented 45 degrees to horizontal and for 
circular polarizations. 5,6 These expectations are supported by mea- 
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surements for terrestrial propagation paths, but there is no experimental 
evidence available to either confirm or refute them for earth-space paths. 
Depolarization caused by ice crystals with unknown orientations raises 
additional uncertainty in the expected behavior of depolarization along 
earth-space paths. 7 

It is difficult and expensive to instrument satellites to measure de- 
polarization directly at several incident polarizations. A technique in- 
volving direct measurement of elements in the polarization transmission 
matrix and direct calculation of depolarization for other incident po- 
larizations was proposed for the COMSTAR beacon propagation experi- 
ments. 8,9 This calculation method has the potential for answering many 
of the questions regarding rain and ice depolarization along earth-space 
propagation paths. This paper presents what are believed to be the first 
experimental results obtained using the technique. 10 A minimum in 
depolarization occurs for vertical and horizontal incident polarizations 
for the attenuation event analyzed. 

II. THE MEASUREMENTS 

2. 1 Equipment and signal parameters 

The transmitting source for the propagation measurements is the 
19-GHz beacon 11 on the COMSTAR satellite located at 95° W longitude. 
This beacon output is switched at a 1-kHz rate between two linear or- 
thogonal polarizations. These 19-GHz signals, among others, are received 
at Crawford Hill, New Jersey with a 7-m diameter antenna and precision 
receiving electronics described in Refs. 12 and 13. 

The polarizations of the beacon signals received at Crawford Hill are 
rotated 21 degrees from horizontal and vertical (for simplicity, these 
signals are referred to as H and V). Amplitudes of copolarized, V and H, 
and cross-polarized, XV and XH, signal components are measured for 
the two transmitted polarizations. Phase differences, </>, referenced to 
the V signal are measured for the H, XV, and XH signals. The receiver 
3-dB predetection bandwidths of 10 Hz in co- and crosspolarized 
channels yield a dynamic range of 60 dB between the clear air copolarized 
signal level and the receiver noise level; crosspolarized signal amplitudes 
are also measured in 1-Hz bandwidths yielding a 70-dB dynamic range. 
Post-detection bandwidths are 1 Hz for all amplitude and phase mea- 
surements. 

Residual signals in the crosspolarized signal channels are produced 
by transmitting and receiving antenna imperfections. These residuals 
are <-35 dB below the copolarized signal levels everywhere within the 
-3 dB beam width of the receiving antenna. Residuals are <— 40 dB 
everywhere within the beam for the XV channel (XV is the channel for 
transmit vertical and receive horizontal); on-axis residuals are typically 
<-45 dB for XV and -36 to -40 dB for XH. 
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2.2 Calibration 

The crosspolarized signal channels are calibrated with the antenna 
feed in the normal receiving position. The receiver polarization switches 
(see Figs. 1 and 14 of Ref. 13) are temporarily configured to connect the 
H and V copolarized signals through the XV and XH channels respec- 
tively. These known signal levels are then used to establish the gains of 
the crosspolarized signal channels. These channel gains remain within 
±0.1 dB over many months. The H and V signals are within 0.05 dB of 
each other as measured by physically rotating the receiving antenna feed 
assembly 90 degrees. The phase scale for the differential phase, 0V-H, 
is calibrated by reversing the polarity of one signal input to produce a 
180-degree phase change. 

Determination of the zero on the phase scales for the crosspolarized 
signals, 0V-XV and $V-XH, requires simultaneous insertion of signals 
with known phase relationship into the receiving feeds. Rotating the feed 
assembly inserts projected spatial components of the copolarized signals 
into both feeds. These components have the same phase and will produce 
signals with either 0- or 180-degree phase difference in the crosspolarized 
signal channels; the sense, or 180, depends on the direction of feed 
rotation. Contamination by the residual crosspolarized signals produces 
an uncertainty in this zero of about ±1 degree. This zero calibration 
remains within ±2 degrees over many months. The phase scale is cali- 
brated by reversing the polarity of one signal input while rotating the 
feed assembly for the zero determination. Note that the phases of the 
clear-air residual crosspolarized signals are not suitable for determining 
phase scale and zero because they are low level and their phases are ex- 
tremely sensitive to the accuracy of the nulling of the residuals with feed 
assembly rotation. 7 

2.3 Baseline removal 

The amplitudes and phases of co- and crosspolarized beacon signals 
exhibit small diurnal variations. These diurnal baseline variations repeat 
very closely from day to day, but change somewhat from season to season. 
The variations result from thermal changes in the satellite beacon cir- 
cuits, waveguides, and antennas as the solar illumination changes. 
Day-to-day receiver contributions to baseline variations are insignificant. 
Because of the repeatability of these baseline variations, they can be 
removed from the experimental data. 

V and H attenuations and differential phase, </>V-H, are determined 
by subtracting, from values measured during propagation events, the 
values measured at the same clock times on clear days within a few days 
of the events. The differential values of attenuation and phase that are 
critical in the calculation of depolarization for other incident polariza- 
tions are determined within ±0.2 dB and ±0.5 degree. 
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Baseline removal of residual crosspolarized components is done by 
vectorially subtracting, from the values measured during events, the 
weighted residual components at the same clock times. These vector 
subtractions can be made since both amplitudes and phases are mea- 
sured. Prelaunch measurements of the beacon antennas, pattern range 
measurements on the receiving antenna, and measurements with the 
receiving feed at normal polarization and rotated 90 degrees all suggest 
that the XH residual is dominated by the beacon antenna. Thus, the XH 
residual is weighted by the event attenuation and differential phase as 
though it were all contributed by the beacon. Weighting of the XV re- 
sidual assumes equal contribution by beacon and receiving antennas. 
The vector subtraction of clear-air crosspolarized channel residuals 
suppresses the residuals to -45 to -50 dB below copolarized signal levels 
before and after events and on adjacent clear days. This accuracy de- 
grades with the intensity of the propagation event because of the un- 
certainty in assigning residual contributions to beacon and receiving 
antennas. However, during weak events when depolarization is low and 
residual suppression is needed, the residual is suppressed to —45 to —50 
dB; when the event is stronger and the suppression degrades, the de- 
polarization is strong and the relative errors contributed by the residual 
are less significant anyway. 

III. POLARIZATION ROTATION TECHNIQUE 

This section introduces the concept of the medium transmission 
matrix for a given set of incident and received polarizations. Receiver 
outputs directly generate the transmission matrix for the incident po- 
larizations; results for other incident polarizations are generated through 
rotations of the measured matrix. 

Let us assume two orthogonal signals, Er a and Erb, incident on the 
medium. The medium outputs in the same reference frame are Er q and 
ERb. The medium is then completely described by four coupling coef- 
ficients, as shown in Fig. 1. Coefficients A and D represent the copolar- 
ized signal attenuations of the medium, while B and C indicate the me- 
dium's depolarization. Note that these coefficients are, in general, 
complex. The above relation may be written in matrix form: 
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Fig. 1 —Description of propagation medium by coupling coefficients A, B, C, D (Era 
and Erb transmitted, Er and E^ received). 
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Since the 19-GHz COMSTAR beacon transmits alternately on two or- 
thogonal polarizations, Er a and Erb are alternately 0. Receiver out- 
puts 8 - 13 are thus directly proportional to A and C (for E T b = 0) or B and 
D (En, = 0). The transmission matrix Tm for the actual incident po- 
larizations is given by these measured values. 

This measured transmission matrix may now be used to generate the 
transmission matrix for any other orthogonal set of incident polariza- 
tions. This operation is done without recourse to path models. Linear, 
circular, or elliptical polarizations are admissible; the following discussion 
considers only linear polarization. As shown in Fig. 2, the transmission 
matrix for the 1-2 frame is desired, based on measurements obtained 
in the M„ — Mb frame. Components in the 1-2 frame are obtained from 
those in the M a — Mb frame through the following relation: 

ca ■ ™ [£]■ 

where 

,_ , fcos 6 —sin 6~\ , , 

[«• * • a A- (4) 

L sin cos v J 

To generate the transmission matrix for incident polarizations along the 
1-2 axes, the 1-2 inputs are rotated into the measurement frame by [Rg], 
passed through the transmission matrix [Tm], and rotated back to the 
1-2 frame by [Re]' l , i.e., 

The transmission matrix [T M ,o] for two orthogonal polarizations rotated 
an angle 6 from the measurement frame is thus 

[TmA = [Re]- 1 [T M ][Ro]. (6) 




Fig. 2— -Relation between reference frames of measured (M a — Mb) and desired (1 - 
2) transmission matrices. 
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IV. OBSERVATIONS 

The accuracy of the polarization rotation technique described here 
is critically dependent on amplitude and phase calibration of all four 
receiver channels. 9 To verify the end-to-end operation of the receiver 
and rotation software, a depolarizer with known symmetry axes was 
inserted near the receiving antenna focus in clear weather. Data were 
taken at several orientations of the depolarizer with respect to the actual 
received polarizations. For each case, depolarization was computed as 
a function of incident polarization angle relative to the depolarizer 
symmetry axes. All results were identical, with minimum depolarization 
along the depolarizer symmetry axes. This agreement validates both the 
receiver calibration and rotation software. 

Data from a recent rainstorm were analyzed to determine the de- 
pendence of depolarization and attenuation on incident polarization 
angle. Since little depolarization was observed before or after the high- 
attenuation portion of the event, it appeared that rain, rather than ice, 
was the predominant depolarizing agent. 7 - 10 

Some effects of incident polarization angle are shown in Fig. 3 for three 
points in this storm. These points had maximum vertical copolarized 
attenuations of 5.4, 10.4, and 20.3 dB. The lower three curves show the 
effect of vertical depolarization (the ratio of vertical crosspolarized to 
copolarized levels) as a function of polarization rotation from local ver- 
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Fig. 3 — Variations in depolarization and differential attenuation with incident polar- 
ization angle for three points during April 20, 1978 rainstorm. Solid curves are calculated 
from measurements made at an incident polarization of -21 degrees indicated by the 
dashed line. 



3272 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1978 



tical to horizontal. The experimental data, taken at a 21 -degree rotation 
angle, appear along the dashed line. All three curves exhibit a sharp null 
within 2 degrees of local vertical, indicating a mean raindrop canting 
angle close to degree. 5,6,10 

The upper three curves indicate the variation in differential attenu- 
ation (the ratio of horizontal to vertical copolarized attenuations) with 
polarization angle. Maximum differential attenuation coincides with 
minimum depolarization, since the raindrops then exhibit their maxi- 
mum oblateness along the two incident polarizations. Differential at- 
tenuation changes sign at 45 degrees from this point, as "vertical" and 
"horizontal" interchange roles. 

The temporal history of a portion of this rainstorm is shown in 
Fig. 4. The upper two curves indicate maximum depolarization (for 
45-degree polarization) and the concurrent copolarized attenuation. A 
maximum depolarization of — 11 dB was observed at 30-dB copolarized 
attenuation. Maximum depolarization exceeded —20 dB for all copo- 
larized attenuations exceeding 11 dB. 

The lower curves indicate time variations for polarization angles in 
the vicinity of the nulls shown in Fig. 3. The central curve indicates the 
polarization angle exhibiting minimum depolarization, while the outer 
curves show the —20 and —30 dB depolarization contours around this 
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Fig. 4 — Temporal history of portion of April 20, 1978 rainstorm. Upper curves show 
copolarized attenuation and maximum depolarization. Lower curves show angular location 
and extent of depolarization minima. 
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minimum. For this portion of this event, the polarization angle exhibiting 
minimum depolarization remained within 2 degrees of local vertical. 
Depolarization remained below - 30 dB for all polarizations within 
3.5 degrees of local vertical, and below -20 dB within 12 degrees of local 
vertical. 

The relation between attenuation and depolarization appeared to be 
well-behaved for the section shown of this rainstorm. Other regions, 
however, appeared less homogeneous. Figure 5 shows the polarization 
angle dependence of depolarization and differential attenuation for two 
points in the storm separated by 19 minutes. Both points had 10.5-dB 
copolarized attenuation. The earlier point, however, had 2.5-dB lower 
maximum depolarization, almost 50 percent lower differential attenu- 
ation, and more than 50 percent greater maximum differential phase 
shift between copolarized channels. In addition, the depolarization null 
was filled in below —41 dB. All these changes are consistent with lower 
average raindrop ellipticity (or less preferential drop orientation) and 
the presence of a small amount of depolarization due to ice crystals. 7 ' 10 
Very preliminary observations of other events have suggested that de- 
polarization nulls occasionally deviate far from local vertical during 
periods of predominantly ice depolarization. 
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Fig. 5 — Variations in depolarization and differential attenuation with incident polar- 
ization angle for two points in April 20, 1978 rainstorm exhibiting identical copolarized 
attenuation but different depolarization characteristics. Solid curves are calculated from 
measurements made at an incident polarization of —21 degrees indicated by the dashed 
line. 
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V. CONCLUSIONS 

A method has been presented, based on transmission measurements 
at one set of incident polarizations, for direct calculation of attenuation 
and depolarization of the transmission medium for any incident polar- 
ization. This technique has been applied to 19-GHz data from the 
Crawford Hill COMSTAR beacon propagation receiver. Test results in- 
dicate that the receiver calibration accuracy is sufficient to provide useful 
information on polarization-dependent effects along earth-space 
propagation paths. 

The rainstorm analyzed exhibited a sharp depolarization minimum 
for incident polarizations within 2 degrees of vertical and horizontal. This 
implies that the mean raindrop symmetry axis was nearly vertical. The 
minimum depolarization was sufficient for communication systems 
employing polarization reuse. 

Similar observations of ice-produced depolarization events suggest 
that the depolarization minima for these events occur occasionally at 
polarizations other than vertical and horizontal. This implies ice particle 
alignment by other than gravity. The simultaneous presence of ice and 
rain with differing symmetry axes could destroy the depolarization nulls 
observed with either separately. Such events, if frequent, could have a 
serious impact on dual-polarization communication systems. The 
technique described here will allow further investigation of this and other 
polarization-dependent propagation phenomena on earth-space 
paths. 
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